Beyond Miransky Scaling 
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We study the scaling behavior of physical observables in strongly-flavored asymptotically free 
gauge theories, such as many-flavor QCD. Such theories approach a quantum critical point when 
the number of fermion flavors is increased. It is well-known that physical observables at this quantum 
critical point exhibit an exponential scaling behavior (Miransky scaling), provided the gauge coupling 
is considered as a constant external parameter. This scaling behavior is modified when the scale 
dependence of the gauge coupling is taken into account. Provided that the gauge coupling approaches 
an IR fixed point, we derive the resulting universal power-law corrections to the exponential scaling 
behavior and show that they are uniquely determined by the IR critical exponent of the gauge 
coupling. To illustrate our findings, we compute the universal corrections in many-flavor QCD with 
the aid of nonperturbative functional renormalization group methods. In this case, we expect the 
power-law scaling to be quantitatively more relevant if the theories are probed, for instance, at 
integer Nt as done in lattice simulations. 



I. INTRODUCTION 



Strongly-flavored asymptotically free theories, such as 
QCD and QED3 are currently very actively researched. 
In particular, QCD with many quark flavors has drawn 
a lot of attention in recent years. On the one hand, the 
number of (massless) fermions can be considered as an 
external parameter. Such gauge theories are then ex- 
pected to exhibit a quantum phase transition from a chi- 
rally broken to a conformal phase when the number of 
fermion flavors is increased. On the other hand, the un- 
derstanding of strongly-flavored gauge theories underlies 
(walking) technicolor-like scenarios for the Higgs sector, 
see e. g. Refs. [2-|i|. 

The phase structure of gauge theories with Nf fermions 
can indeed be rich, as simple considerations may already 
suggest. Due to the screening property of fcrmionic fluc- 
tuations, asymptotic freedom is lost for large Nf. For 
instance, SU(7V C ) gauge theory with Nf fermions is no 
longer asymptotically free (a.f.) for Nf > N* := ^-N c - 
Another special fermion number Np BZ potentially ex- 
ists denoting the minimum flavor number for the occur- 
rence of an infrared fixed point g 2 of the running gauge 
coupling. For instance, the two-loop /? function of the 
gauge coupling g 2 exhibits the so-called Caswell-Banks- 
Zaks (CBZ) fixed point as the screening nature 

of fermion fluctuations dominates the two-loop coeeffi- 
cient for Nf > Np BZ . For instance for SU(3), we have 
Nf Z ~ 8.05 in the two-loop approximation. A pertur- 
bative treatment of the theory seems possible near A r f a f -, 
Nf < iV f •, where g 2 is small, indicating the existence of a 
conformally invariant limit in the deep infrared [ll[ . For 
decreasing Nf, g 2 becomes larger, suggesting the onset 
of chiral symmetry breaking. The decoupling of massive 
fermions then destabilizes the Caswell-Banks-Zaks fixed 
point g 2 in the gauge sector of the theory. The infrared 
of the theory is then dominated by massless bosonic ex- 
citations, the Goldstonc modes, and the spectrum of 



the theory is characterized by a dynamically generated 
mass gap . A similar reasoning also applies to QED3, see 

c g. Gllll. 

These considerations suggest the existence of a quan- 
tum critical point associated with a critical flavor num- 



ber A f CBZ < A fcr < 



jya.f. 



above which gauge theories 



approach a conformally invariant limit in the infrared. 
Thus, Nf serves as a control parameter for the quantum 
phase transition. 

Studies of the phase structure of strongly-flavored 
gauge theories have been performed employing contin- 
uum methods as well as lattice simulations. In QED3 
many studies have performed estimates to determine 
JVf jCr using Dyson- Schwinger equations and resumma- 
tion techniques [r3 - l22j . Since the dynamically generated 
mass is substantially smaller than the scale set by the 
gauge coupling, lattice simulations of QED3 with many 
flavors are remarkably challenging |23T - |26| . The phase 
structure of many-flavor QCD has also been studied em- 
ploying continuum methods El, I27l - l43| . as well as 
lattice simulations [H|-[58| . Recent results suggest in this 
case that a conformal phase indeed exists with a quantum 
phase transition occurring near 9 < A f cr < 13. 

Given the existence of such a quantum critical point 
in an asymptotically free gauge theory with Nf flavors, 
the question arises how the spectrum of the theory be- 
haves when we approach this quantum critical point 
from below. This question is tightly bound to the ques- 
tion of the Nf dependence of the dynamically generated 
scale associated with chiral symmetry breaking. It is 
well-known from studies of Dyson-Schwinger equations 
in the rainbow-ladder approximation that physical ob- 
servables, e. g. the fermion condensate, exhibit an ex- 
ponential scaling close to Nf tCI , provided that the (mo- 
mentum) scale dependence of the gauge coupling can be 
neglected HIH-Ell, 



fc SB oc A0(iVf, cr - Nf) exp - 



2ey/ ai \N{ tCC -N f \ 



(1) 
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Here, ksB denotes a scale characteristic for the onset of 
symmetry breaking, being directly proportional, say, to 
a symmetry-breaking condensate. The quantities e and 
<X\ are pure constants arising from the details of the 
theory and will be defined below. This behavior can 
be viewed generalization of essential Berezinskii- 

Kosterlitz-Thouless (BKT) scaling [62M64} to higher di- 
mensional systems [65(. We rush to add that the spec- 
trum of the different theories below and above iVf iCr are 
substantially different. In particular, a construction of an 
effective low-energy theory in terms of light scalar fields 
may no longer be possible above iVf jCr . 

Taking into account the running of the gauge coupling 
and going beyond the standard rainbow-ladder approxi- 
mation, the scaling behavior of physical observables close 
to iVf.cr has been analyzed in [33, [H, [66| . More precisely, 
the Nf dependence of a strict upper bound for the sym- 
metry breaking scale has been studied, scaling according 
to a power law, 

fc cr oc A\N fiCr -N{\~^. (2) 

Here, fc cr denotes the scale where the dynamics leading 
to symmetry breaking becomes critical. This means that 
operators that trigger symmetry breaking become rele- 
vant in an RG sense. As the system still has to run 
towards lower energy scales into the broken phase, we 
have k cr > fcsB , implying that Eq. ^ is an upper bound 
for Eq. (UJ. Near the critical flavor number the corre- 
sponding scaling exponent is uniqely determined by the 
critical exponent Oo of the gauge coupling at its infrared 
fixed point. This upper bound for the (chiral) sym- 
metry breaking scale can be translated into an upper 
bound for physical observables (6(|. In fact, the chiral- 
phase-transition temperature as a function of the "exter- 
nal" control parameter Nf has been computed with non- 
pcrturbative functional renormalization group methods. 
The scaling of the phase boundary has been found to be 
compatible with the analytically derived scaling behav- 
ior jMm. 

Recently, the scaling behavior of physical observables 
has been investigated again with the aid of Dyson- 
Schwinger equations in the rainbow-ladder approxima- 
tion also taking into account part of (momentum) scale 
dependence of the gauge coupling by a proper adjustment 
of the scale [67| ■ It was then found that the exponential 
scaling behavior close to iVf jCr of the symmetry breaking 
scale is supplemented by a power-law behavior similar to 
the one found in Rcf . [37], [H, HH ■ 

In the present work, we aim to reveal the relation be- 
tween these supposedly different scaling laws and show 
rigorously what kind of scaling behavior we should ex- 
pect close to the quantum critical point of asymptotically 
free gauge theories with many flavors. Our arguments 
are based on very general RG considerations and involve 
only a few assumptions about the underlying fixed-point 
structure of the theory under consideration. In fact, we 
shall show that the above-given scaling laws arise as two 
different limits of one and the same RG flow. In addition, 



we point out the importance of the scale-fixing procedure 
applied in the first place in order to compare theories 
with different flavor numbers. As the scaling behavior 
of the low-energy observables is accessible to a variety of 
nonperturbative methods, we believe that a rigorous un- 
derstanding of scaling behavior near the phase transition 
to the conformal phase will be very useful. 

In Sect. UH we briefly repeat the arguments that lead 
to an exponential scaling behavior at the quantum phase 
transition. In addition, we derive the leading-order cor- 
rection to the exponential scaling behavior. In Sect. IIII1 
we then discuss the issue of scale fixing underlying a 
meaningful comparison between theories with different 
flavor numbers. Moreover, we briefly review the argu- 
ments from Refs [13, HE which lead to a power-law- 
like scaling behavior for a strict upper bound for the 
(chiral) symmetry breaking scale. In Sect. IIVI we then 
discuss the interrelation of the scaling laws put forward 
in Refs. HE HE HE H3 an d derive the leading-order scal- 
ing behavior of a given infrared observable at the quan- 
tum critical point. To illustrate our analytic findings, we 
present numerical results from a non-perturbative func- 
tional renormalization group study of the scaling behav- 
ior in many-flavor QCD in Sect. [Vj 



II. MIRANSKY SCALING 

In this section we study exponential scaling behavior in 
gauge theories near a quantum critical point, also known 
as Miransky scaling [23, H9| . 

We shall keep our discussion as general as possible. 
For our purposes, however, we consider a general class of 
theories where symmetry breaking and condensate for- 
mation is driven by fermionic self-interactions. Indepen- 
dently of whether these interactions may be fluctuation- 
induced (as in QCD) or fundamental, this class of theo- 
ries can be parameterized by the following action: 

Sm = J d d x^ip(i$ + g4)i> + Kp-rS^aippipyips} ,(3) 

where a, (3 ... denote a specific set of collective indices 
including, e. g., flavor and/or color indices. In general, 
we expect to have more than just one four-fermion in- 
teraction channel, see e. g. Sect. [V] for QCD with many 
flavors. Note that symmetry breaking is ultimately trig- 
gered by the interactions approaching criticality, i.e., be- 
coming RG relevant. 

Here and in the following we do not allow for terms 
in the action which explicitly break the underlying 
symmetry, such as current quark mass terms in QCD 1 . 



The scaling behavior of observables with the (current) quark 
mass in the (quasi-)conformal phase of strongly-flavored gauge 
theories is of particular interest for lattice simulations and cur- 
rently under investigation, see Refs. [68147011 . 





(a) (b) (c) 

Figure 1: Representation of the terms on the right-hand side of the RG flow equation Q by means of Feynman diagrams. Our 
functional RG studies, see Sect.|Vl include resummations of all diagram types including ladder-diagrams generated by type (b) 
and (c) as well as the corresponding crossed-ladder topologies. 



From the action © we can derive the /3 function of 
the dimensionless four-fermion coupling A. It assumes 
the following simple form: 



Px = d t X = (d - 2)A - aX 2 - b\g 2 



eg 



(4) 



Here, t = ln(fc/A) denotes the RG 'time' with k be- 
ing the RG scale and A being a UV cutoff scale. The 
couplings A ~ A/fc( d-2 ) and g ~ g/k 4 ~ d denote dimen- 
sionless and suitably rcnormalizcd couplings. The first 
term in Eq. ([4]) arises from simple dimensional reseat- 
ing. The quantities a, b and c do not depend on the 
RG scale but may depend on control parameters, such as 
the number of quark flavors N{ or the number of colors 
A^. 2 The various terms on the right-hand side of Eq. (|4]) 
can be understood in terms of perturbative Feynman di- 
agrams [7l|, see Fig.Q] Equation ((H) can also be derived 
from nonperturbative flow equations in the limit of point- 
like (momentum-independent) interactions, sec Sect. [V] 
Note that we have dropped terms proportional to the 
anomalous dimension of the fermionic fields in Eq. 
We assume these contributions to be small in the follow- 
ing. This is indeed the case in the chirally symmetric 



72] 



regime of QCD, at least in the Landau gauge 

In Eq. © we have not further specified the gauge sec- 
tor. In fact, let us ignore the running of the gauge cou- 
pling in this section, and consider the gauge coupling as 
a scale-independent "external" parameter. The RG flow 
of the gauge coupling is then trivially governed by 



d t g 2 = . 



(5) 



For example, this may be an acceptable approximation in 
the vicinity of an IR fixed point g 2 . Note that the value 
of g 2 may depend on other control parameters such as 



2 Note that the coefficients a, b and c can depend implicitly on 
the RG scale as soon as we introduce a dimensionful external 
parameter, e. g., temperature T. However, the coefficients re- 
main dimensionless since they depend only on the ratio T/fc, see 
e. g. MM- 



Nf or AT C , cf. our discussion of QCD with many flavors 
in Sect. 

In Fig. [5] we show a sketch for the j3\ function, implic- 
itly assuming that a > 0, b > and c > in Eq. (QJ. For 
a vanishing gauge coupling g 2 we find two fixed points, 
an IR attractive Gaussian fixed point at A = and an 
IR repulsive fixed point at A > 0. For increasing g 2 these 
fixed points approach each other and eventually merge 
for a critical value g 2 r , 



d-2 



6 + 2 



(6) 



For g 2 > g 2 r the four-fermion coupling then becomes 
a relevant operator and increases rapidly towards the 
IR indicating the onset of (chiral) symmetry breaking. 
Thus, the four-fermion coupling A necessarily 3 diverges 
for g 2 > g 2 r at a finite RG scale kss = fcsB(<? 2 )- This 
divergence is, of course, an artifact of the over-simplistic 
approximation ([3]), but can be related to a symmetry- 
breaking transition in the effective Landau-Ginzburg- 
type potential for fermion-bound states. Even though 
ks~B is not a direct observable, it sets the scale for ob- 
servables such as condensates, decay constants, critical 
temperatures, etc. This picture of the emer gen ce of chi- 
ral symmetry has been put forward in (3(14381 [66j and 
successfully employed for an anlysis of the phase struc- 
ture of QCD with various numbers of flavors and colors 
at zero and finite temperature (36M381 [66j . Moreover, 
this picture has also been employed to study conformal 
scaling in quantum field theories, see e. g. Ref. jr35j . 



3 Here, we assume that the initial conditions at the UV scale k = A 
for the four-fermion coupling A are chosen such that Aa is smaller 
than the value of the IR repulsive fixed point, see Fig. [2] In 
beyond-standard model applications Aa is sometimes considered 
to be a finite parameter, see e. g. [73l . We therefore add that 
the exponential scaling behavior discussed below can only be 
observed when Aa is chosen to be smaller than the value of the 
repulsive fixed point for a given q 2 . Otherwise, we expect a 
power-law-like scaling behavior l74t . 



Let us now briefly discuss the scaling behavior of the 
symmetry-breaking scale k$B when g 2 is varied by hand 
as a constant "external" parameter. To this end, we have 
to solve the RG flow equation (QJ. We find: 



2 arctan 



In k - In A = - - 



bg 2 -(d-2)+2aA' 
<5(9 2 ) 



% 2 ) 



with 



S(g 2 ) = ^Aacg 4 -{{d-2)-bg 2 ) 2 . 



(7) 



(8) 



Here and in the following we assume b > without loss 
of generality. From Eq. ([7]) , we obtain fcsB by solving for 
the zero of 1/X(k), i. e. 1/A(fcge) = 0: 



In k 



SB 



In A 



S(9 2 ) 



const. . 



(9) 



Here, we have chosen the initial conditions such that 



A 



uv 



A max where A r 



denotes the position of the max- 



imum of the P\ function, i.e., the peak of the parabola 
in Fig. [3J An expansion of © around g 2 T yields 4 



fc SB oc A6(g 2 - g 2 r ) exp 



where e is a numerical factor, 



2^(g 2 -g 2 CI ) ' 



(10) 



(d- 2) (2ac + by/ac) 
b + 2^/ac 



(11) 



which in general depends on the details of the theory un- 
der consideration, e. g. the number of colors and flavors 
in QCD, In any case, we find an exponential Miransky- 
scaling behavior of ksB for g 2 close to g 2 r . Since the 
dynamically generated scale fcsB sets the scale for the 
low-energy sector, we expect that all IR observables O 
scale according to: 



= fc 



^SB 1 



(12) 



where do is the canonical mass dimension of the observ- 
able O and fo is a function which does not depend on 
g 2 T but may depend on g 2 and other external param- 
eters, e. g., Nf and/or N c . The function fo can be 
computed systematically within certain approximations 
schemes such as large- N c expansions or chiral perturba- 
tion theory, see e. g. (6(| H3]- In the context of QCD 
the scaling law in Eq. (|TU|) has been first derived by Mi- 
ransky 27, H^] but has also been found in the context 
of specific 2-dimensional condensed-matter systems [HJ . 



Note that g 2 r is defined to be the value of g for which the f3\ 
function has exactly one zero. In general there exist two solutions 
for g 2 T , however, one of which can be excluded from a physical 
point of view. 




Figure 2: Sketch of a ty pical /3 function for the fermionic self- 
interactions A, see [3(| and also [38[ for a generalization to 
finite temperature): at zero gauge coupling, g 2 — (upper 
black curve), the Gaufiian fixed point A = is IR attractive. 
For g 2 = g 2 r (middle/blue curve), the fixed-points merge due 
to a shift of the parabola induced by the gauge-field fluc- 
tuations ~ g 4 . For gauge couplings larger than the critical 
coupling g 2 > g 2 r (lower /red curve), no fixed points remain 
and the self-interactions rapidly grow large, signaling chiral 
symmetry breaking. The arrows indicate the direction of the 
flow towards the infrared. 



The derivation of the scaling law (jTU)) via an analysis of 
the RG flow of four-fermion operators has been recently 
pointed out by Kaplan, Lee, Son and Stcphanov [65 1. 

Let us now briefly discuss the consequences of the 
scaling law (|10l) when we apply our considerations to 
strongly-flavored gauge theories, such as QCD with many 
quark flavors or QED 3 . In these cases we may choose the 
IR fixed-point of the gauge coupling as an external pa- 
rameter, i. e. g 2 = g 2 {Nf) in Eq. (fTUj) . Depending on 
the Nf dependence of the coefficients a, b and c in the /3a 
function, the critical value for the gauge coupling may 
depend on the number of flavors as well, g 2 r = g^^Nf). 
The critical number of quark flavors Nf lCT can then be 
obtained from the criticality condition 



ffc 2 r (iV fiCr ) = «72(iV f , cr ). 



(13) 



This corresponds to the coupling value for which the 
two fixed points of the four-fermion coupling A merge 
and then annihilate each other for g 2 > g 2 T . Expanding 
gl{Nf) - gl{Nf tCI ) around N {<CI , 

9 2 *(Nf)-gi.(Nf. cr ) = a 1 {N i -Nf. CI ) + a 2 (Nf-N lcl .) 2 + .. . . 

(14) 

and plugging (fl"4"l) into (fTU)) . we find the exponential Nf 
scaling of k$B'- 



k SB oc A9(Nf tCr -N { ) exp 



7T(1 



\N f . 



-N { \ + 



2e^/ ai \N lcT -N { \ 



(15) 

Whether the size of the regime for exponential scaling is 
small depends on the ratio a^/cti which in turn depends 
on the theory under consideration. Thus, the size of the 
scaling regime may presumably be different in, c. g., QCD 
and QED 3 . In Sect. [V] we compare the analytic findings 
of this section with results from a numerical analysis of 
QCD with many flavors. 
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III. POWER-LAW SCALING 

In this section wc discuss how the running of the gauge 
coupling affects the RG flow of the four-fermion cou- 
pling^). In particular, we argue that (chiral) symmetry 
breaking in strongly-flavored gauge theories is a multi- 
scale problem, in contrast to the scenario discussed in 
the previous section. In other words, the (chiral) sym- 
metry breaking scale fcsB discussed above and its scaling 
with the control parameters, e. g. the number of flavors 
Nf, depends on the scale fixing and its potential flavor 
dependence. 

In the following, we include the running of the gauge 
coupling which goes beyond standard rainbow-ladder ap- 
proaches employed in the context of strongly-flavored 
gauge theories, see e. g. Ref. (32j . 

Since we arc eventually interested in the scaling behav- 
ior of IR observables, e. g. the fermion condensate, it is 
important to realize that a variation of the flavor number 
does not quite correspond to a change of a parameter of 
the theory. It rather corresponds to changing the theory 
itself. We would like to stress that there is indeed no 
unique way to unambiguously compare theories of differ- 
ent flavor number with each other, as different theories 
may be fixed at different scales. 

As we have argued in more detail in Ref. |66j, fixing 
the scale of theories with, say, different flavor numbers Nf 
by keeping the running coupling at some scale A (e. g. 
r mass) fixed to a certain value, seems to be a well ac- 
cessible prescription for many non-perturbativc methods. 
In general, it is important to take care that this scale- 
fixing procedure is not (or as little as possible) spoilt by 
scheme dependencies. The latter constraint essentially 
rules out Aqcd as a proper scale in QCD to be kept 
fixed in theories with different flavor numbers. Of course, 
it is also possible, c. g. in lattice QCD simulations, to 
keep the value of an IR observable fixed for theories with 
different Nf, e. g. the pion decay constant or the crit- 
ical temperature. We shall briefly comment on such a 
procedure below. For what follows, however, we choose 
a mid-momentum scale for the scale fixing, lying in be- 
tween the high-scale perturbative running and the more 
interesting non-perturbative dynamics. Thus, wc fix the 
theories at any Nf by keeping the running coupling at 
some intermediate scale A fixed to a certain value, say 
a \. 

To be more specific, wc shall focus our discussion on 
strongly-flavored asymptotically free gauge theories, such 
as QCD with many flavors and QED 3 . 5 In such theories, 
the dependence of the running coupling on the scale and 
on further control parameters such as Nf is expected to 
modify Miransky scaling. In particular, an understand- 
ing of the universal scaling behavior of observables in the 



a By asymptotic freedom, we refer to the vanishing of the dimen- 
sionless renormalized coupling in the UV. 



ordered phase close to the phase transition at iVf iCr is 
of interest. However, the arguments also apply to other 
theories in which dynamical chiral symmetry breaking is 
trigged by a running coupling which approaches a non- 
trivial IR fixed point. 

For a monotonically increasing coupling flow, the value 
of the non-trivial IR fixed point gl of the gauge coupling 
corresponds to the maximum possible coupling strength 
of the system in the conformal window, i. e. for Nf >cr < 
Nf < Nf A \ As both gl and g 2 T depend on the number 
of flavors, the criticality condition gl(Nf iCr ) = g 2 r (Nf yCT ) 
defines the lower end of the conformal window and thus 
the critical flavor number, see Sect. Inland the left panel 
of Fig. [3] for an illustration. 

For gl > g 2 r , our model ^ is below the conformal 
window and runs into the broken phase. Slightly below 
the conformal window, the running coupling g 2 exceeds 
the critical value while it is in the attractive domain of 
the IR fixed point g~. The flow in this fixed-point regime 
can approximately be described by a (3 function expanded 
around the fixed point g 2 : 

/3 g2 ee d t g 2 = -Q (g 2 - g 2 ) + 0((g 2 -gl) 2 ) . (16) 

The universal "critical exponent'" denotes (minus) the 
first expansion coefficient. We know that 6 < 0, since the 
fixed point is IR attractive, see right panel of Fig. [3] In 
general, the critical exponent depends on Nf, 9 = Q(Nf). 
The solution to Eq. (fTBf for the running coupling in the 
fixed-point regime simply reads 

3 2 (fc)^-(£)"°, (17) 

where the scale k is implicitly defined by a suitable ini- 
tial condition and is kept fixed in the following as we keep 
the UV scale A fixed. 

The scale fco corresponds to a scale where the sys- 
tem is already in the fixed-point regime. For the present 
fixed-point considerations, fco provides for all dimension- 
ful scales. However, from the knowledge of the full RG 
trajectory, ko can be related to the initial scale A, say 
the r mass scale in QCD, by RG evolution. 

A necessary condition for (chiral) symmetry breaking 
is that gl > g 2 r . This implies that g 2 (k) exceeds g 2 T at 
some scale fc cr which is implicitly defined by the criticality 
condition, g 2 (Nf tCV ) — g 2 r (N{ yCV ), and therefore 

fccr > fc'SB , (18) 

where fcse is the scale at which the four-fermion coupling 
A diverges, see Sect. [H] Thus, fc cr is an upper bound for 
the symmetry breaking scale fcsB ■ From Eq. (|17|) and the 
criticality condition g 2 (k CI ) = g 2 r , we derive an estimate 
for fc cr valid in the fixed-point regime 

k CI ~k (gl~g 2 cr )-^. (19) 

The scale fc cr is dynamically generated. Note that 
fccr/fco — > for gl —> g 2 r from above. Due to our scale- 
fixing procedure, this scale depends on Nf and Nf iC1 - in a 
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Figure 3: Left panel: illustration of the IR running of the gauge coupling in comparison with the critical value of the gauge 
coupling g 2 r . Below the conformal window, Nf < JVf. cr , g 2 exceeds the critical value g 2 r , triggering the approach to xSB. For 
increasing flavor number, the IR fixed-point value g 2 becomes smaller than the critical value, indicating that the theory is 
inside the conformal window. Right panel: sketch of the /3 function of the running gauge coupling. The slope of the j3 function 
at the IR fixed-point corresponds to minus the critical exponent 0, cf. Eq. (|16l) . The vertical line to the right gives the value 
of g 2 r . The dotted line gives the value of the gauge coupling at the UV scale A which we keep fixed for all Nf. By contrast, the 
value of g 2 T may depend on Nf. The arrows indicate the direction of the flow towards the infrared. 



non-trivial way 6 . Using Eq. (|14[) and a Taylor expansion 
of the critical exponent near the quantum phase transi- 
tion, 

0(iV f ) = e + 0i(JV f - N f , a ) + <D((Nf - AT fiCr ) 2 ) , (20) 
we find the following Nf dependence of k CI for Nf < Nf a - 

k c , ~ fcb|JV fl0r -JVfr^r (21) 

/ \N [iCr -Nf\ fa 2 e x \\ 

+ ... , 

where 0o = 0(7Vf jCr ). Since k cr defines the scale at which 
the fixed-points in the j3 function of the four-fermion cou- 
pling merge, the existence of a finite k cr can be considered 
as a necessary condition for (chiral) symmetry breaking. 
Thus, we expect that the scale for a given IR observables 
O for Nf < iVf )Cr is set by fc cr : 

O = fokif , (22) 

where do is the canonical mass dimension and fa depen- 
des on Nf but not on Nf iCI , see also Eq. (fl2|). However, 
we would like to stress that k CT does not include the full 
dependence of fcge on (Nf — iVf jCr ), i. e. k cr /ksB ^ const, 
is still a function of the control parameter, as we shall 
discuss in the subsequent section. 

In Ref. [33, [H, [6(| we have implicitly used the existence 
of a finite k cr to estimate the chiral phase transition tem- 
perature in QCD as a function of Nf. For a given value 



of iVf the phase transition temperature is given by the 
highest temperature for which we still have fc cr > 0. We 
have indeed found that T cr scales according to Eq. (j21j) : 



T cr ~ k \N fiCI -Nf\ 4. (23) 



Strictly speaking, this is only an upper bound for the chi- 
ral phase transition temperature since it is only sensitive 
to the emergence of a fermion condensate on interme- 
diate (momentum) scales but insensitive to a fate of the 
condensate in the deep IR close to T cr due to fluctuations 
of Goldstone modes [75(, also known as a local ordering 
phenomena. Such strong IR fluctuations of the Gold- 
stone modes may yield corrections to the scaling law for 
the critical temperature given above 7 . Nevertheless, re- 
lation (|23|) is an analytic prediction for the shape (of the 
upper bound of) the chiral phase boundary in the (T, iVf) 
plane. 

At vanishing temperature, the analysis of the scaling 
behavior of IR observables is simplified compared to a 
scaling analysis at finite temperature since dimensional 
reduction does not set in in the deep IR enhancing the 
Goldstone modes. Based on the observed scaling behav- 
ior of fc cr with the number of flavors, we are therefore in 
a position to derive the Nf scaling of low-energy observ- 
ables, such as fermion condensates, at zero temperature. 



6 Note that it is, in principle, possible to adjust the initial value 
of the coupling at the initial scale such that the scale k CI is 
independent of Nf and Nf CT . As indicated above, we expect 
that such a scale-fixing procedure would, however, be strongly 
affected by scheme-dependencies at least in our truncation. 



7 We would naively expect that corrections to Eq. I|23|l can be 
only resolved in lattice simulations with very small masses for 
the pseudo Goldstone modes and on very large lattice sizes, see 
also Sect. |IV] 
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IV. BEYOND MIRANSKY SCALING 

Let us now discuss how the symmetry breaking scale 
^sb < ^cr depends on (N { — Nf cr ). We consider again a 
Lagrangian of the form and assume that Nf < Nf CI . 
The crucial new ingredient compared to the derivation 
of Miransky scaling is the RG flow of the coupling. We 
also assume that the system has already evolved from 
the initial UV scale A to the scale fc cr at which the fixed 
points of the f3 function of the four-fermion coupling have 
merged. Sufficiently close to -/Vf lCr , the flow of the gauge 
coupling is governed by the fixed point regime for g 2 > 
g 2 T . The running of the gauge coupling is then given by 
(cf. Eq. (H3) 



9 2 (k) = gl-{ 9 l- 9 l) 



9l - (Aff 2 ) 



k 



(24) 



where Ag 2 = g 2 — g 2 x . Recall that g\ ~ Nf and Ag 2 
\Nf tCI - Nf\. Plugging Eq. into Eq. (0J, we find 

-e 

p x = d t \ = P x (Ag 2 



ij- 



(d-2)X-aX 2 - bXg 



eg* 



k 

n cr 

dg 2 



(25) 



k 



where we have used Eq. (|T9|) . Recall that k < k cr ^ ko 
and O < 0. We observe that the zeroth order in Ag 2 
coincides with the f3\ function for which we have found 
an (implicit) analytic solution for constant g 2 in Sect. [TT] 
yielding Miransky scaling. We refer to this analytic solu- 
tion as \ g 2 . The solution of the (3 function (f25j) can then 
be found by an expansion around the solution X g 2 : 



X = X g 2 + (Ag 2 )SX + . . 



sx 



(26) 



This allows us to systematically compute the scaling be- 
havior for Nf < -/Vf. cr . Since we are interested in the 
(chiral) symmetry breaking scale k^B we have to solve 
1/A(fcge) = for fcgB- In zeroth order the scale fcge can 
be computed along the lines of our analysis in Sect. [TXJ 
We find 

fc SB oc k CI 6(Nf iCI - N f ) cxp -- — 1 = ] 

V 2ey'ai|iVf ;Cr -7Vf|y 

~ k 6(Nf, CI - Nf)\Nf, CI - Nf\~^ 
/ 7T \ 



x exp 



2e^/ai|iVf, cr -iVf|y ' 



(27) 



where we have used Eq. (|2ip in leading order. Higher 
order corrections to Eq. (|27[) can be computed system- 
atically as outlined above and in the previous sections. 



Thus, we have found a universal correction to the expo- 
nential scaling behavior which is uniquely determined by 
the universal "critical" exponent 9. A similar result has 
been suggested very recently by Jarvinen and Sannino us- 
ing a standard rainbow-ladder approach with a constant 
gauge coupling but a properly adjusted scale [67[. Our 



RG analysis demonstrates in a simple and systematic way 
that such a rainbow-ladder approach is indeed justified 
and yields the correct leading-order scaling behavior. 

Let us now turn to the scaling behavior of physical ob- 
servables. The scale of all low-energy observables is set 
by &sb- In other words, k^n represents the UV cutoff of 
an effective theory at low energies, such as chiral pertur- 
bation theory, quark-meson or NJL-type models in case 
of QCD. At zero temperature we therefore expect that 
a given IR observable O with mass dimension do scales 
according to 



O = fo{Nf)6{Nf <cx - Nf) 



(28) 



ffM) - VNf 



where /o(iVf) is a function which depends on Nf but 
not on Nf_ cr . As mentioned above, fo(N) can be in 
principle systematically computed in QCD using, e. g., 
chiral perturbation theory or a large-A^ c expansion. For 
instance, in a large- iV c expansion it is straightforward to 
derive the leading Nf dependence of the function fo (Nf) 
for the pion decay constant /„.. It reads (66| 



(29) 



The scaling law (|28j) together with (|27)) represents one 
of the main results of this work. It can be used as an 
ansatz to fit, e. g., data from lattice simulations. This 
scaling law is remarkable for a number of reasons: first, it 
relates two universal quantities with each other: quanti- 
tative values of observables and the IR critical exponent. 
Second, it establishes a quantitative connection between 
the (chiral) phase structure and the IR gauge dynamics 
(0). Third, it is a parameter-free prediction following 
essentially from scaling arguments. Moreover, it shows 
that Miransky scaling and power-law scaling are simply 
two limits of the very same set of RG flows: in the limit 
— > oo we find pure Miransky-scaling behavior, while 
we have pure power- law scaling in the limit 0^0. 

At this point, we would like to emphasize once more 
that the scaling behavior of any IR observable near iVf cr 
depends crucially on the scale-fixing procedure applied 
in the first place. Still, the universal scaling will always 
show up at one or the other place and thus cannot be re- 
moved, as stressed in Ref. [66| ■ Our choice to fix the scale 
at m T which is large enough not to be affected by chiral- 
symmetry-brcaking is certainly not unique. In principle, 
the point where to fix the scale can be chosen as a free 
function of Nf. In Eq. (TIT)) , this would correspond to 
the choice of an arbitrary function ko = ko(Nf) for the 
global scale, which then appears also in the scaling re- 
lations (|2"Tj) . and (|2~7} . Indeed, an extreme choice 
would be given by measuring all dimensionful scales in 
units of a scale induced by chiral symmetry breaking 
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(such as T cr or /„■). In this case, all chiral observablcs 
would jump non- analytically across Nf = N". Neverthe- 
less, the scaling relations would then translate into scal- 
ing relations for other non-chiral external scales: e.g., the 
scale k at which the running coupling acquires a specific 

value (say a = 0.322) would diverge with TVf — > iV f cr ac- 

_ i 

cording to k ~ \Nf — iV f cr | e oi . This point of view can 
constitute a different way of verifying our scaling rela- 
tions on the lattice. 

Let us conclude this section with a discussion of the 
importance of the corrections to the exponential scaling 
behavior due to the running of the gauge coupling. To 
that end, it is convenient to consider the logarithm of the 
(chiral) symmetry breaking scale /csb, 

1 7T 

In fc SB = const. - — In |JV f , cr - N { \ - - — ===== . 

6 2ey/ ai \N ljCI -N f \ 

(30) 

This expression can be used to estimate the regime in 
which the corrections to the exponential scaling become 
subdominant. For this, we compute the minimum of the 
function 



1 



I6n 



In \Nf, a - N f \ + 



(31) 



with respect to |iVf iCr — iVf |. In accordance with Eq. (f2"U)) . 
we assume \Nf iCt — Nf\ < 1 here. From this, we can 
then estimate that corrections to the exponential scaling 
behavior are subdominant as long as 



|iV f -iV fjCr | 



< 



^ 2 |©o| 2 
16e 2 ai ' 



(32) 



with e begin defined in Eq. (fTTj) . We observe that cor- 
rections to Miransky scaling due to the running of the 
gauge coupling are small when |0o| 3> 1 and large when 
1 0o | <C 1. To be more specific, let us consider QCD with 
many flavors: assuming N{ >CT ~ 12, we extract 0o ~ 0.4 
from the two-loop f3 g 2 function. From Eq. (|32[) the region 
where pure Miransky scaling dominates is then found to 
be confined to the regime |iVf — iVf lC r| ^ 0.3. Thus, we ex- 
pect that the exponential scaling behavior is dominantly 
visible only very close to -/Vf jCr . According to this esti- 
mate, the ©-dependent universal corrections are there- 
fore more significant in QCD. 

In QCD, it appears to be a general feature that ©o 
decreases with iVf. cr . Estimates of ©o within two- and 
higher-loop approximations in the MS scheme are sum- 
marized in Fig. |4] Therefore, power-law scaling is more 
prominent for larger Nf tCI . In particular, power-law scal- 
ing should be visible if theories are probed only for integer 
values of Nf as, e.g., on the lattice. 

The role of |0| for the scaling behavior close to -ZVf jCr 
can also be understood by simply looking at the f3 g 2 func- 
tion of the gauge coupling, see Fig. [3l For |0| 3> 1 the 
gauge coupling runs very fast into its IR fixed point once 
it has passed g 2 x . Thus, the situation for g 2 > g 2 r is as 
close as possible to the situation studied in Sect. |TTJ The 




Figure 4: Critical exponent <d of the running gauge coupling 
at the Caswell-Banks-Zaks fixed point as a function of the 
number of flavors Nf as obtained from two-, three- and four- 
loop perturbation theory in the MS scheme. 



coupling can simply be approximated by a constant. For 
|0| <C 1 the gauge coupling runs very slowly ("walks") 
into its IR fixed point once it has passed g 2 T . This walking 
behavior for g 2 > g 2 r then gives rise to sizable corrections 
to the exponential scaling behavior. 

Finally, we would like to discuss the finite-temperature 
many- flavor phase boundary in QCD. In [37l l38j it was 
found that the scaling of the phase boundary is consistent 
with the pure power-law scaling behavior (|23p . From the 
above discussion this result is now understandable since 
the exponential scaling behavior sets in only very close to 
-/Vf.cr for iVf cr w 12 and thus remains invisible in numer- 
ical fits over a wider range of Nf. Of course, power-law 
scaling behavior for the chiral phase-transition temper- 
ature still remains an upper bound even if we took into 
account the exponential factor in Eq. (|28|) . This is due to 
the fact that strong fluctuations of Goldstone modes in 
the IR may yield further corrections and lower the phase 
transition temperature, see e. g. Ref. [75j . Whether these 
corrections at finite temperature yield additional correc- 
tions to the scaling behavior cannot be answered within 
the scaling analysis presented in this work. However, it 
may very well be that such corrections depend only on 
Nf but not on iVf jCr . Nevertheless, we would like to stress 
that a further investigation of the finite-temperature scal- 
ing behavior at the quantum critical point, Nf = Nf fCr , 
seems worthwhile in QCD since the scaling behavior in 
Nf direction may significantly differ from the expected 
power-law scaling behavior in the temperature direction 
at fixed Nf, see I76{ . 



V. QUANTITATIVE SCALING ANALYSIS IN 
MANY-FLAVOR QCD 

Having derived analytic scaling relations for physical 
observablcs in the previous sections, we present results 
from a numerical study of the scaling behavior in QCD 
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with many flavors in this section. 



A. Renormalization group setup 

Our numerical analysis is based on previous works 
on strongly-flavored gauge theories in the framework of 
a functional RG approach using the Wetterich equa- 
tion [77], for reviews, see 



In [36[ the zero- 
temperature quantum phase transition of QCD with iV c 
colors and Nf flavors has been computed using the func- 
tional RG. The phase diagram at finite temperature as 
a function of Nf has first been computed in [37l HI] . We 
briefly review these results in this section and employ 
them for our numerical analysis of the scaling behavior. 

In [36l - [38j the RG flow of QCD starting from the micro- 
scopic degrees of freedom in terms of quarks and gluons 
was studied within a covariant derivative expansion. A 
crucial ingredient for chiral symmetry breaking are the 



scale-dependent gluon-induced quark self-interactions of 
the type included in © . We note that dynamical quarks 
influence the RG flow of QCD by qualitatively different 
mechanisms. First, quark fluctuations directly modify 
the running of the gauge coupling due to the screening 
nature of these fluctuations. On the other hand, gluon 
exchange between quarks induces quark self-interactions 
which can become relevant operators in the IR as we 
have already discussed in the previous sections. These 
two mechanisms strongly influence each other as well. 
As we have seen, however, it is possible to disentangle 
the system once we accept that these fluctuations can be 
associated with different scales in the problem. 

From now on we restrict ourselves to d = 4 Euclidean 
space-time dimensions and work solely in the Landau 
gauge. In a consistent and systematic operator expan- 
sion of the effective action, the lowest nontrivial order is 
given by [72j 



J 



d i xl$(tf + g4)il> + - A_(V-A) + A + (V+A) + A CT (S-P) + A VA [2(V-A) ad j+(l/iV c )(V-A)] \. (33) 



r 



This ansatz for the effective action underlies our non- 
perturbative RG study. The four-fermion interactions oc- 
curring here have been classified according to their color 
and flavor structure. Color and flavor singlets are 



(V-A) 
(V+A) 



(V'7/xV') 2 



(^7m75^) 2 



(34) 
(35) 



where (fundamental) color . . . ) and flavor £, . . . ) 
indices are contracted pairwise, e.g., (ipip) = {tj)fipf). 
The remaining operators have non-singlet color or flavor 
structure, 



(S-P) : 
(V"A) adj : 



(^)M^W) 2 ^(#^) 2 - 
(^t>) 2 + (^ 7m75 t>) 2 , 



(V?75^) 2 
(36) 



where ('tp x ipt) 2 = 2p x tp^-i/j^tp x , etc., and {T a )ij denote 
the generators of the gauge group in the fundamental 
representation. 

We stress that the set of fermionic self-interactions 
introduced in Eq. (|33|) forms a complete basis. This 
means that any other pointlike four-fermion interaction 
which is invariant under SU(iV c ) gauge symmetry and 
SU(7Vf)L x SU(7Vf)R flavor symmetry can be related to 
those in Eq. ([33]) by means of Fierz transformations. In 
our numerical analysis, we neglect UA(l)-violating inter- 



actions induced by topologically non-trivial gauge con- 
figurations, since we expect them to become relevant 
only inside the x^B regime or for small Nf. In addi- 
tion, the lowest-order UA(l)-violating term schematically 
is ~ {^>iji) Nt . Thus, larger Nf correspond to larger RG 
"irrelevance" by naive power-counting. Moreover, inter- 
actions of the type ~ ('ifjifj) 1 * 1 for Nf > 3 do not contribute 
directly to the flow of the four-fermion interactions due 
to the one-loop structure of the underlying RG equation 
for the effective action. 

As a severe approximation, we drop any nontrivial mo- 
mentum dependencies of the A's and study these cou- 
plings in the point-like limit A(|p;| <C k) in our scaling 
analysis. Therefore our ansatz for the effective action 
does not allow us to study QCD properties in the chirally 
broken regime, since, e. g., mesons manifest themselves 
as momentum singularities in the A's. Nonetheless, our 
point-like approximation can be reasonable in the chi rally 
symmetric regime. This has been indeed shown in [36J, 
where the regularization-scheme independence of univer- 
sal quantities has been found to hold remarkably well in 
the point-like limit. 



Using the truncated effective action (|3"3"|) . we obtain 
the following j3 functions for the dimensionless couplings 
A, = Ai/fc 2 (see |sSlz3): 
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d t X- = 2A- 



a ;( FB )> 4 



A 7 , 



■,g 2 A_ - 3# 2 Ava 



1 ,(FB),4 



12 + 9iV e 2 

7V2 ' 



-%v A lf )A { - NfN c (\ 2 _ + A 2 ,) + A 2 _- 2(/V c + iV f )A_A V A + A^A+A, + 2A 2 /A }, 



(37) 



d t X+ = 2A + -4w 4 4 F 1 B) ' 4 

-8^ F) < 4 {-3A 2 + 



3 2, 



1 ,(FB),4 



12 + 37V 2 



(38) 



2JV c JV f A_A + - 2A+(A_ + (N c + N { )X VA ) + N { X-X a + X VA X a + iA ff 2 }, 



d t X„ 



2X CT -Av i lZi [6C 2 (N c )g 2 X a - 6g 2 X + ] - -v 4 l\ 



(FB),4 
2 



^8 U4 4 F),4 {2iVcA 2 r -2A_A a -2^ f A ff AvA-6A + A CT }, 



24 - 9N 2 
N 



(39) 



ft A 



'(AVA 



2AvA-4 t ;44 F 1 B) ^ 4 



N 



g 2 X V A - 3g 2 A. 



1 ,(FB),4 
o V * l l,2 



24 - 37V 2 

Wr. 



-8v, 



4 1 { P A [-(N C + 7V f )A 2 /A +4A_A VA - W 2 }. 



r 



(40) 



Here, C 2 (N C ) = (N 2 - 1)/(2N C ) is a Casimir operator of 
the gauge group, and V4 = l/(327r 2 ). The regularization- 
schcmc dependence of the RG flow equations is con- 
trolled by (dimcnsionless) threshold functions I which 
arise from Feynman diagrams and incorporate fermionic 
and /or bosonic fields |Sa|. For the optimized regula- 
tor I89h91|. we find 



(F),4 



1 ,(FB),4 _ 1 VA ,(FB),4 _ 3 

2' 11 ~ 6 ' h ' 2 ~ 2 



VA 

6 



(41) 



In our numerical analysis, we have dropped contribu- 
tions from the anomalous dimensions of the fermions and 
the gauge coupling t]a = (3 g 2/g 2 . While the first one is 
proportional to the gauge-fixing parameter and vanishes 
identical in the Landau gauge in the chirally symmet- 
ric regime [721 ] . we have found by a comparison of our 
numerical results with those from [36j that the contri- 
butions oc rjA in the threshold function do not strongly 
affect our result for Nf >CT , In fact, we have t\a —> for 
N f 



iV f a f - and g 2 < g 2 . Moreover, we find for g 2 < g 2 



< 1 for Nf > 11 and \v A ' l °° P \ ^0-5 for 
In total, this may lead to quantitative correc- 



that 

Nf > 

tions at most on the percent level. 

Let us now discuss the running of the gauge coupling. 
Even though the running coupling has been computed 
within the functional RG approach [3?], HU, [92ll94 |. we 
employ for simplicity the two- and four-loop result ob- 
tained in the MS scheme [H, ISfiJ , as our results show a 
satisfactory convergence in the strongly flavored regime. 
We will often restrict ourselves to the two-loop case, as 
it already shows all qualitative features and can be dealt 



with analytically. The analytic expression for the two- 
loop /3 g 2 function reads explicitly: 



with 



d t g 2 



11 



( 9 2 



V16tt 2 



9 

8tt 2 



Po = ^rN c - -N h A = 



34iV c 3 + 3N { - 13N 2 N { 



(42) 



(43) 



Note that the chosen regularization scheme in the matter 
sector and the MS scheme do not coincide. This inconsis- 
tency results in an error for our estimate for the critical 
number of quark flavors. Since we are rather interested in 
the scaling behavior which is related to the universal crit- 
ical exponent 0, our results are only influenced indirectly 
by this approximation 8 . Due to this scheme dependence, 
the results using the four-loop running may not necessar- 
ily be considered as a more precise calculation. Instead, 
the difference between two-loop and four-loop MS results 
should be viewed as an estimate of the dependence of our 
results on the quantitative details of the running gauge 
sector. 

A comment on contributions to the running of the 
gauge coupling induced by the presence of the quark 
self-interactions Ai is in order here: To render the RG 



8 Of course, the actual value of ©0 = ©(-Wfcr) depends on the 
actual value of JVf C r which itself, as a universal quantity, depends 
on the difference of the scheme-dependent quantities g^ r and gl . 
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theory. For our numerical study we choose the fixed- 
point value of the gauge-coupling at the two-loop level: 



2 

9*, 2- 



Am 



16(lliV c 2 - 27V c 7vf)7r 2 



(44) 



-loopV^ -i/ l 3N 2 N{ _ UN 3 _ 3 JV f 

In the matter sector we employ two different trunca- 
tions to which we refer as one-channel and all-channels 
approximation. The latter one is Fierz complete. In the 
all-channels approximation we take into account the full 
set of flow equations (f5T]) - (|4T)]) . while we only take into 
account the RG flow of the scalar-pseudoscalar channel 
A (j in the one-channel approximation and set all other 
four-fcrmion couplings to zero: 



Figure 5: Logarithm of the (chiral) symmetry-breaking scale 
\n(ksB/m T ) as a function of (iVf jCr — N{ tCI )/N{ tCI for an Nf- 
dependent but scale-independent, i.e. constant gauge cou- 
pling. The corresponding fits are given in Eq. (f52)) . 



flow gauge invariant we have to take regulator-dependent 
Ward-Takahashi identities into account [H, H?} . In the 
present case, these symmetry constraints yield contribu- 
tions to the running of the gauge coupling which de- 
pend on the quark self-interactions. However, these 
contributions are proportional to the j3 functions of 
the four-fermion couplings, as has been pointed out in 
Ref. [36|, [72|. Therefore, these contributions vanish as 
long as the four-fermion couplings are at their fixed 
points, i. e. as long as g 2 < g 2 r . Thus, we expect that 
these contributions do not alter the scaling-law ([27)) in 
leading order 9 . In particular, the power-law behavior is 
unaffected by these corrections arising due to symme- 
try constraints. In the present approximation, we ignore 
these corrections in our numerical analysis. 



B. Miransky-type scaling 

Let us start with a numerical analysis of many-flavor 
QCD with a constant gauge coupling: 

d t g 2 = . 

As discussed above, the gauge coupling can then be con- 
sidered as an "external" iVf-dependent parameter of the 



cag 



(45) 



with 



Nc_ , _3_ 

4tt 2 ' 4tt 2 
3 /97V 2 -24 



C 2 {N C ) . 



256tt 2 



N, 



(46) 



Here, we have adopted the conventions of Sect. [IT] for the 
coefficients a, b, c. 

Now we can compute the critical values of the gauge 
coupling in the one- and in the all-channels approxima- 
tion. In the one-channel approximation we find 



2 

"cr.one 



32tt 2 (W c 3 - 27V c - ^3iV c 6 - 87V 4 ) 



(N c =3) 



10.E 



3(4 + iV 4 ) 

(47) 

which does not depend on Nf. In the all-channels ap- 
proximation the critical value has to be computed nu- 
merically. As found in (36|, the resulting critical value 
9cr all 01 t ne gauge coupling then depends on Nf, for a 
given number of colors, g 2 r all decreases weakly with in- 
creasing Nf. 

The fixed-point value g 2 2 _i oop together with the crit- 
ical value of the gauge coupling can be used to estimate 
the critical number of quark flavors above which there 
is no chiral symmetry breaking in the IR. In agreement 
with the results given in Ref. [36j , we find 



J 



N? nc 

lv f.cr 



169iV c 6 - 136iV 4 + 132A 2 - 68^ N* (3iV 2 -8)iV c 3 
58iV c 5 - 647V 3 - 26 ^A 4 (3A 2 -8)A 2 + 6^A 4 (37V 2 -8) + 36 N c 



(N c =3) 



11.7, 



(48) 



r 



for the one-channel approximation and 
N?» w 11.9 



(49) 



for the all-channels approximation. We may use our es- 
timate for JVf jCr from the one-channel approximation to 
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estimate N{ cr in the limit A c — >• oo: 



N onc 
f,cr 



68^- 169 
2(13^3-29) 



3.95. 



(50) 



Our results for Af )Cr arc in accordance with the results 
from Dyson-Schwinger equations in the rainbow-ladder 
approximation, see e. g. [a, IH, 01 as well as with those 
from current lattice simulations [4J-|57[ . 

Let us now study the dependence of the symmetry 
breaking scale fcge on Af for the specific case A c = 3. 
In Fig. [S]we show our results for ln(fcgB/A) as function 
of (Af jCr — Nf) /Nf as obtained from the one-channel (dots) 
and from the all-channels (triangles) approximation us- 
ing gl 2-ioop as a fixed input parameter. As initial con- 
ditions for the Aj's for a given g\ 2 -ioo P (-^f) we have used 
the solution of the coupled set of linear equations 



dXi 



= 0. 



(51) 



where i £ {+,— ,er, VA}. This corresponds to starting 
the flow at the maxima (extrema) of the parabolas. 

We observe that for a given Nf the symmetry break- 
ing scale fcsB is smaller in the all-channels approximation 
compared to the one-channel approximation. The fits to 
the data points are also shown in Fig. [5] In agreement 
with our analytic results we find: 



In k 



sl; 



const. — 



2.481 



|A f>cr - iV f | - 494 ' 
Inkfi « const.- |Ar 3 ' 9 ^, n , B1fi • ( 52 ) 



iiv f ,cr - A^r 



Thus, we clearly observe the expected exponential scal- 
ing behavior in the one-cannel and in the all-channels 
approximation for Nf Af jCr . 

The result from the one-channel approximation is in 
reasonable agreement with the analytic leading-order 
(LO) result found in Sect, fll] 



InfcLO 



const. 



const. 



2e^K||A f , cr 

2.386 
y/\N f , CI - A f | 



A f | 



(53) 



Note that |a2/ai| « 0.273. Differences to the numerical 
results are due to numerical errors of the fit and higher- 
order corrections which we have derived in Sect. [TTJ 



C. Power-law scaling and beyond 

Let us now take into account the (momentum) scale- 
dependence of the running gauge coupling. In order to 
compare the theories with different flavor numbers we fix 
the scales by keeping the running coupling at the T-mass 
scale A = m T fixed to a(m T ) « 0.322. Since we ap- 
ply the truncation (|3"3")l to QCD, we do not consider the 
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Figure 6: Nf dependence of fc cr and ksB as obtained from 
a study with a running coupling in the two-loop approxima- 
tion. The criticality scale fc cr (blue circles) is dominated by 
power-law scaling (straight line with slope ~ | ©o | 1 in this 
double-log plot), and clearly serves as an upper bound for the 
symmetry breaking scale fcsB (red triangles), being a super- 
position of power-law and Miransky scaling. If the theories 
are probed at integer Nf, i.e., AiVf > 0(1), the contribution 
due to Miransky scaling may not be visibile. A pure powerlaw 
fit to chiral observables ~ &sb , may however overestimate the 
critical exponent Oq. 



four-fermion couplings A as independent external param- 
eters as, e.g., in Nambu-Jona-Lasinio-type models. More 
precisely, we impose the boundary condition Xi — > for 
k — > oo which guarantees that the A's at finite k are solely 
generated by quark-gluon dynamics, e.g., by 1PI "box" 
diagrams with 2-gluon exchange, cf. Fig. [He). 

In Fig. [6] and Fig. [7] we show our results for the Nf 
dependence of the scales fc cr and ksB as obtained from 
a study with a running gauge coupling in the two-loop 
and the four-loop approximation, respectively. Note that 
Af. cr becomes smaller when we employ the running cou- 
pling in the four-loop approximation. We find 



N 



4 — loop 



10.0 



(54) 



4— loop 
f,cr 



in the all-channels approximation and N t 
in the one-channel approximation, in agreement with 
Ref. [H. 

The data points can be fitted to our analytic results for 
the scaling behavior of fege and k cr . For the all-channels 
approximation, we find 

\nk^- loop w const. + 2.566 |JV f -JV f , cr |, (55) 

3.401 



and 



In k 2 ~ loop 
m k sb 



lnfci- loop 



const. — 



In k 



4— loop 
SB 



\Nf- A f , cr |0-54 
-2.540 In I Af- A f , cr | , 



const. + 1.180 | A f - A fjCr | , 
5.196 

const. 



|A f - A f , cr |0-52 
1.171 In I A f - A f , cr |. 



(56) 
(57) 

(58) 
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Figure 7: N{ dependence of fc cr and fcsB as obtained from 
a study with a running gauge coupling in the four-loop ap- 
proximation. The contributions due to Miransky scaling, 
roughly parameterized by the difference between k CI (blue 
circles) and fcsB (red triangles), extend to larger values of 
ANf = 7Vf iCr — Nf, as the estimate for the critical exponent 
Bo = 0(A r f jCr ) at four loop is larger than at two-loop. In this 
perturbative estimate for the running coupling, the curves 
cannot be extended to larger values of AJVf, see text. 



Thus, the fits are in reasonable agreement with our an- 
alytic predictions. For the multi-parameter fits (|56[) 
and ([55]). we have fixed the coefficient of the ln-term 
which is the inverse critical exponent On = Q(Nf CI ). It 
should be stressed that the predicted values for the crit- 
ical exponent Q(Nf. cr ) are substantially different for the 
running coupling in the two- and four-loop approxima- 
tion (cf. Fig 13} , 



1 



|e(jv f , cr )| 
i 

|0(^cr)| 



2.540 (two- loop) , 

1.171 (four-loop). (59) 



In Fig. IS] and Fig. [7] we observe that the critical expo- 
nent O clearly influences the scaling behavior close to 
the quantum critical point Nf^ cr . In agreement with our 
analytic findings, the size of the regime with exponential 
scaling increases with increasing critical exponent 0. Us- 
ing Eq. (|32[) we can give a quantitative estimate for the 
size of the regime in which the exponential scaling be- 
havior dominates. For the one-channel approximation, 
cf. Eq. (gSJl, we find: 

ANf := \Nf - iV f , cr | < 0.3 (two-loop, Nf, cr 11.7). (60) 

Using a running coupling in four-loop order (Nf lCT 9.8) 
the size of this Miransky scaling regime can be estimated 
to be larger than one flavor. This is in agreement with 
our numerical results, see Figs. [5] and [7] With this per- 
turbative estimate for the running coupling, however, the 
curves in Figs. [6] and [7J cannot be extended to larger val- 
ues of ANf = N{ >CI — Nf. For instance, in the four-loop 
case, we have Nf >CI ~ 9.8 on the one hand. On the other 
hand, the Caswell-Banks-Zaks fixed point vanishes for 



Nf < 8. Our RG arguments based on expansions about 
an IR fixed point hence only extend to AiV"f >max ~ 1.8, 
cf. Fig. [7] In nonperturbative functional studies where 
an IR fixed point appears to exist already in the pur e 
gauge sector and thus also at lower Nf [33, HH, I98l4l05 1 , 
no restriction on ANf arises. 

To summarize: since iVf iCr > 9 in current lattice simu- 
lations (44T - [57| , we expect that the pure exponential scal- 
ing behavior is difficult to resolve and the corrections 
due to the running of the gauge coupling (0) might be 
more relevant for lattice simulations. From the viewpoint 
of such simulations, one might be interested in keep- 
ing the power of the "Miransky" term fixed to 1/2 and 
use the scaling law to fit Nf cr and the critical exponent 
0o = B(Nf iCI ). 

To illustrate the influence of the critical exponent we 
have also computed the scaling behavior of the scales fc cr 
and ksB using a model for the running gauge coupling. 
This model is inspired by the two-loop approximation 
modified by an artificial higher-order term. The latter 
is constructed such that the critical exponent can be 
changed by hand, still leaving the two-loop fixed point 
unaffected: 



-loop 



5*,2- 



loop/ 



(61) 



where the parameter <fi allows us to change without 
changing A/f jCr . In Fig. [5] we present our results fcsB and 
k cr for (j) = 6.003 (i. e. |9(iV f , cr )| « 4.3) in the left panel 
and for = -0.0001 (i. e. |6(iVf )Cr )| 0.3) in the right 
panel. The results clearly confirm that the size of the 
exponential-scaling regime depends strongly on 0. 



VI. CONCLUSIONS 

In this work we have analyzed how physical observ- 
ables in asymptotically free gauge theories, such as QCD 
or QED3, scale when the number of flavors is varied. 
When the number Nf of fermion flavors in such theories 
is increased, a regime may open up along the Nf axis in 
which the theory is asymptotically free but remains chi- 
rally symmetric in the infrared. This gives rise to the 
existence of a quantum critical point on the Nf axis. The 
exact determination of the location of this quantum crit- 
ical point in QCD as well as in QED3 is currently a very 
active frontier in theoretical physics. 

Even though we have presented estimates for Nf rr as 
obtained from a functional RG approach, see also [3a . l3a [. 
the focus of this work is on the actual scaling behavior of 
physical obscrvables close to the quantum critical point. 
This scaling behavior of observables such as the fermion 
condensate close to -/Vf jCr is not only interesting in its own 
right but may also help to guide future lattice simulations 
in this field. 

Ignoring the running of the gauge coupling, it has 
long been known that physical observables obey Miran- 
sky scaling, i.e., an exponential scaling law close to the 
quantum critical point [27l l59l |62T - [64| . In more recent 
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Figure 8: Nf dependence of k cr (blue circles) and fcsB (red triangles) as obtained from a study with a model for the running 
gauge coupling, cf. Eq. (|61[) . that allows to tune the critical exponent Q by hand. We show the results for So = |0(iVf lCr )| ~ 4.3 
(left panel) and 0o| ~ 0.3 (right panel). Contributions due to Miransky scaling are visible as deviations from a straight-line 
behavior (power law) in this double-log plot. These results confirm our estimate that the Miransky-scaling window is larger for 
larger |0o| (left panel), whereas power-law scaling dominates for small |Go| (right panel). 



studies [33, |66| we have shown that an upper bound for 
the scaling behavior of, e. g. the chiral phase transi- 
tion temperature, can be derived from an analysis of the 
fixed-point structure in the matter sector. In combina- 
tion with the running of the gauge coupling in its fixed- 
point regime, the upper bound for physical observablcs 
then scales according to a power law. The associated 
critical exponent is related to the universal critical expo- 
nent O at the IR fixed-point of the gauge coupling, see 
Eq. (EH). 

From another viewpoint, we have shown that Miransky 
scaling close to iVf iCr receives universal power-law correc- 
tions, which are uniquely determined by the critical ex- 
ponent at the IR fixed-point of the gauge coupling. Both 
scaling laws follow from one and the same set of RG flow 
equations and can be considered as two different limits 
of a general scaling law: pure exponential Miransky scal- 
ing arises in the limit of large O — > 00, whereas power- 
law scaling becomes more prominent at small — > 0. 
Quantitatively, we have estimated the size of the regime 
with almost pure exponential scaling in strongly-flavored 
QCD and found it to be small, |AT f)Cr - N { \ < 0.3 for 
-^Vf,cr ~ 11-7. Outside this regime the scaling behavior of 
physical observablcs is controlled by the critical exponent 
0. Our numerical analysis of scaling in many-flavor QCD 
based on functional RG methods is indeed in agreement 
with these analytic findings. 

Finally we would like to add that the scaling behavior 
close to the quantum critical point can be contaminated 
by the scale-fixing procedure. As also argued in [66|, a 
comparison of theories with different Nf is not unique 



for non-conformal theories but indeed requires a specific 
choice of a dimensionful scale. This scale is used as a ruler 
for the different theories. In the present work, we have 
fixed the scales by keeping the gauge coupling fixed to the 
same value for all Nf at an initial mid-momentum scale, 
c. g. the r mass scale in QCD. Alternatively, theories 
with different Nf can be fixed by keeping an IR observable 
characteristic for the ordered phase for all Nf fixed, say 
the pion decay constant f n or T CI in QCD. However, the 
behavior of physical observables at the quantum critical 
point is then discontinuous. 

In any case, the scaling relation (|28[) is a parameter- 
free testable prediction for the behavior of physical ob- 
servables near the quantum critical point. On the one 
hand, it might be tested directly by lattice simulations. 
On the other hand, our prediction might also be helpful 
to guide future lattice simulations of strongly-flavored 
gauge theories. 
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